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Abstract. We study the tangential case in 2-dimcnsional almost- Riemannian 
geometry. We analyse the connection with the Martinet case in sub-Rieman- 
nian geometry. We compute estimations of the exponential map which allow 
us to describe the conjugate locus and the cut locus at a tangency point. We 
prove that this last one generically accumulates at the tangency point as an 
asymmetric cusp whose branches are separated by the singular set. 



1. Introduction 

In a series of recent papers [2,3,8], 2-dimensional almost- Riemannian geometry 
is investigated under generic conditions, giving rise to Gauss-Bonnet type results 
on compact oriented surfaces. 

Roughly speaking, an almost-Riemannian structure (ARS for short) on an n- 
dimensional manifold M can be defined locally by the data of n vector fields playing 
the role of an orthonormal basis. Where the vector fields are linearly independent, 
they define a Riemannian metric. But the structure is richer along the set Z where 
they are linearly dependent (see section [3] for a precise definition of ARS). 

For 2-dimensional ARS, it was proven in [2] that generically the singular set 
Z is an embedded submanifold of dimension I and only 3 types of points exist: 
the ordinary points where the metric is Riemannian, the Grushin points where the 
distribution A generated by the vector fields has dimension I and is transverse to 
Z, and the tangency points where A has dimension I and is tangent to Z. 

The situation around ordinary and Grushin points is well known from the metric 
point of view, even if new considerations about curvature close to the Grushin points 
allow the authors to prove new results in [2,3,8]. 

These metrics have also been studied in [5]. In that paper, the authors deduce 
a global model on the two-sphere of revolution S 2 as a deformation of the round 
sphere, the metric being 

gx = d<p 2 + G x {tp)d6 2 , Ae [0,1], 

with G\(X) = jz^jx, where X — sin 2 tp, and (tp,9) are the spherical coordinates. 
In this representation, the singularity is located at the equator: <p> = n/2. This 
metric appears in orbital transfer and, moreover, the homotopy is important to 
understand the behavior of the curvature. In this framework a short analysis tells 
us that for the generic model the symmetries (of revolution and with respect to the 
equator) cannot be preserved and a non integrable model is obtained. 
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In this paper we analyse the situation at tangency points. The presence of these 
points is fundamental in the study of 2-dimensional ARS. 

In [3] , the authors provide a classification of oriented ARS on compact oriented 
surfaces in terms of the Euler number of the vector bundle corresponding to the 
structure (see[5]for definition) in presence of tangency points, generalizing a result of 
[2] . The construction of Gauss Bonnet type formulae is more intricated in presence 
of tangency points because of the geometry of the tubular neighborhoods of Z close 
to the tangency points (see [3]). 

It happens that the geometry close to tangency points is not well known and 
more intricated for many reasons. First, the computation of expansions of the wave 
front is more complicated and involves elliptic functions. Second, the nilpotent 
approximation is far from being generic as defined below. In particular the dis- 
tribution of the nilpotent approximation is transversal to its singular set at the 
tangency point. 

In this paper we focus on two points. First, we analyse the connection be- 
tween tangency points in 2-dimensional ARS and Martinet points in 3-dimensional 
sub-Riemannian structures. This allows us to obtain regularity properties of the 
distance function. Second, we compute the jets of the exponential map which al- 
lows to estimate the conjugate locus and the cut locus at the tangency point. In 
particular we prove that, differing from the nilpotent case, the cut locus generi- 
cally accumulates at the tangency point as an asymmetric cusp whose branches are 
locally separated by the singular set Z (see figure [TJ) . 




FIGURE 1 . The sphere (solid line) and the cut locus (dashed line) at a 
tangency point in the generic case together with the singular set (dotted 
line) 

The paper is organised as follows. In section [5] we recall some basic definitions 
and results. In section[3]we show the relation between ARS and constant rank sub- 
Riemannian structures. In section 0] we analyse the special case of the nilpotent 
approximation as well as a generic model at a tangency point. In section [5] we 
compute the asymptotic expansions of the exponential map at a tangency point. 
This allows us to estimate, in section [BJ the conjugate and cut loci at a tangency 
point, giving rise to a geometric interpretation of the first invariants in terms of the 
form of the cut locus. 

2. Basic definitions 

An n-dimensional ARS is the data of a triple (M , E, /) where M is an n- 
dimensional manifold, E is a Euclidean bundle of rank n over M and / is a morphism 
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of vector bundles between E and TM preserving the basis M, such that the eval- 
uation at any point q <E M of the Lie algebra generated by {/ o <j \ a section of E} 
is T q M. 

From the control theory point of view, an ARS can be defined locally by the data 
of n vector fields (Fi,..., F n ) such that Lie{Fi, . . . , F n } q = T q M for all q. They 
define locally the following control dynamical system 

n n 

(21) $ = 5>? = i> 

i=l i=l 

the distance between two points qa and qi being by definition the minimal time 
needed to join qi from q n with this control system. We also define the submodule 
A of the module of vector fields on M generated locally by (F\, . . . , F n ) and the 
flag A fc by Ax = A, A fe+1 = A fe + [A, A fc ]. 

In the following we deal with 2-dimensional ARS. Let us recall the following 
result proved in [2]. 

Proposition 1. The following properties, denoted by (H0) 7 are generic for 
2-dimensional ARSs. 

(1) The singular set Z is a one- dimensional embedded submanifold of M , 

(2) the points q € M where A 2 (q) is one- dimensional are isolated, 

(3) A 3 (<7) = T q M for all q G M. 

Moreover, if a 2-dimensional ARS satisfies (HO), then for every point q € M there 
exist a neighborhood U of q and an orthonormal frame (F\,F 2 ) of the ARS on U 
such that, up to a change of coordinates, q = (0,0) and (Fi,F 2 ) has one of the 
forms 

(Fl) F 1 {x,y) = &, F 2 (x,y) = e*^£, 

(F2) F 1 (x,y) = £, F 2 (x,y) = xe^£, 

(F3) F 1 (x,y) = £, F 2 {x,y) = (y-x^{x))e^y)^ 

where <f), ip and £ are smooth functions such that 0(0, y) = and tfj(0) > 0. 

Remark. In order to get the same notations as in [7], the normal form (F3) will be 
written in the following in coordinates (y, z) 

(F3) F 1 (y,z) = %, F 2 {y,z) = {z-y^{y))e^f z . 

For a 2-dimensional ARS satisfying (HO), we say that a point q is 

• ordinary if A(q) = T q M (normal form (Fl)), 

• a Grushin point if the dimension of A(q) is one and A 2 (g) = T q M (normal 
form (F2)), 

• a tangency point if the dimension of A(q) = A 2 (q) is one and A 3 (q) = T q M 
(normal form (F3)). 

In the normal form (F3), (y,z) is a privileged coordinate system with weights 
respectively 1 and 3 (for definitions of privileged coordinates and nilpotent approx- 
imation we refer the reader to [4]). 

Consider the change of coordinates y = y and z = —z/(2ip(0)e^). Let us 
still denote (y, z) by (y, z). According to the weights the jet up to order of the 
elements of the orthonormal frame in the normal form (F3) is 

(F3) F 1 (y,z) = ^ j , F 2 (y,z) = (ez+£+e'y 3 + o 3 (y,z))£ z 

where e — e^ ^ 0, s' = ^ ^^2^(0) dv ^ an< ^ °3(J/j z ) IS a smooth function of order 
higher than 3 in the variables y and z with respect to their weights. 
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3. ALMOST-RlEMANNIAN GEOMETRY AND SUB-RlEMANNIAN GEOMETRY 

3.1. Local desingularization of an n-dimensional ARS. Let us present a clas- 
sical construction. Consider an n-dimensional ARS on M and let {F\, . . . , F n ) be a 
local orthonormal frame on a neighborhood of q. Assume that Fi(q) = for i > d 
where d = dim A(q). Define 

M = M x R n - d = {{x,y)\x e M,y e R"~ d }. 

Denote by 7Ti and 7T2 the canonical projections on the first and second factor of M. 
Then we can define Ff by 

7ru(§)=Ji 

7r 2 *(F 2 ) = if * < d, 

Then the family {.Fi, . . . , i 7 ^} has rank n in the neighborhood of q and defines the 
orthonormal frame of a sub-Ricmannian metric on M. Moreover, if A is bracket 
generating as a submodule of the Lie algebra of vector fields on M, then the same 
holds true for A = span{Fi, . . . , F n }. This metric on M is invariant with respect 
to translations in R n ~ d . Moreover, one can show that the curves between q and 
qi minimizing the almost-Riemannian distance on M are projections of the curves 
between {qo} x H n ~ d and {qi} x R™~ d minimizing the sub-Riemannian distance 
on M, a curve in M and its projection having the same length. This implies that 
the ball centered at q of radius r in M is the projection of any sphere of radius r 
centered at a point of the type (q, y) in M. Applying the Pontryagin Maximum 
Principle to the corresponding extremals in M, the transversality conditions to 
{g } x R"~ d and to {qx} x R"~ d must be satisfied. 

3.2. Examples. 

3.2.1. The Grushin plane and the Heisenberg group. The 2-dimensional ARS de- 
fined by the orthonormal frame |Fi = J| , F 2 — xJ^ j in R 2 is the Grushin plane. It 

is the first example of almost-Riemannian structure with non empty singular locus, 
namely the y-axis. Moreover, it is the nilpotent approximation at any Grushin 
point of a 2-dimensional ARS (for a precise definition of nilpotent approxima- 
tion of a system we refer the reader to [4]). If we apply the desingularization 
procedure, we find the sub-Riemannian metric defined by the orthonormal frame 
{-^i = Jx' ^ 2 = x ~§y + ~§z \ on R- 3 j which is the Heisenberg metric. The Heisen- 
berg metric is the nilpotent approximation at any point of contact of a rank-2 sub- 
Riemannian structure defined on a 3-dimensional manifold, that is at any point p 
where the rank-2 distribution satisfies [A, A](p) = T p M. The Hamiltonian associ- 
ated with the Grushin metric is 

Hi = \(pl + x 2 p 2 y ) 
while the one related to the Heisenberg metric is 

H2 - ^(pl + (xp y + p z ) 2 ) 

where p x , p y , p z are the dual coordinates to x, y and z in the cotangent bundle. 
Geodesies of the Heisenberg group projecting to geodesies of Grushin are those 
with p z = 0, respecting the transversality condition to the vertical lines given by 
the Pontryagin Maximum Principle. 

In general, the relation between the cut and conjugate loci for the sub-Rie- 
mannian metric on M and almost-Riemannian one on M is not clear, the projection 
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tti introducing singularities. The Grushin plane is a good illustration of this fact. 
The two geodesies for the Grushin metric starting at (0, 0) with the initial covectors 
(p x = l,p y ) and (p x = -l,p y ) are 

*(t) = P*^4 
Py 

_ 2p y t - sm(2p y t) 

Hence, these two geodesies first intersect for p y t — ir and one can prove that i 
corresponds to the cut time along them. Moreover, computing the Jacobian of the 
exponential mapping, one proves that the conjugate time t satisfies p y t = tan(p y t). 
Lifting to the corresponding geodesies in the Heisenberg space starting at (0, 0, 0) 
with initial covector (p Xl p y ,p z = 0), one finds as third coordinate 



z(t) = p x 



1 — COS(Pyt) 



Py 

Hence, the two lifted geodesies do not intersect anymore at t and the computation 
of the Jacobian of the exponential mapping shows that the conjugate time for both 
lifted curves satisfies = tan(^). It corresponds to the second conjugate time 
in the Heisenberg case. 

3.2.2. The nilpotent approximation at a tangency point and the Martinet flat case. 
Consider the normal form (F3) at a tangency point as presented in section [2j 
Recall that the weight of the variable y is 1 and the weight of z is 3. Hence the 
nilpotent approximation can be given in (y, z) coordinates by the orthonormal frame 

{ ~§y ' ^Sz }' corresponding to the metric g = dy 2 + (^-)^ 2 dz 2 and the Hamiltonian 



-(p 2 + y -- 2 -- 
2 {Py 4 



Applying the desingularization procedure, one finds the orthonormal frame 

d y 2 d d 
dy : 2 dz dx 

in (y,z, x) coordinates in R 3 , whose corresponding Hamiltonian is 



This lifted structure is the Martinet flat case of sub-Riemannian geometry. It is 
the nilpotent approximation at any Martinet point of a rank-2 sub-Riemannian 
structure defined on a 3-dimensional manifold, that is at any point p where the 
rank-2 distribution satisfies [A, A](p) = A(p) and [[A, A], A](p) = T p M. The set 
of these points is called Martinet surface. Being the nilpotent approximation, the 
Martinet flat case will provide the starting point to analyse the general tangential 
case and it will allow to make some preliminary estimates of the sphere and the 
distance function using previous computations as in [7]. 

4. Local analysis at a tangency point 

In this section we focus on the following models in order to study the local 
situation around tangency points for a generic 2-dimensional ARS. 
(1) The nilpotent approximation (of order -1) 

g _ 1= dy 2 + (^y 2 dz 2 
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(2) The generic model of order 

v 2 

g = dy 2 + (ez + y — + e'y^dz 2 
where e = e' = gives the nilpotent approximation. 

4.1. Analysis of the nilpotent model. In this case, the desingularization pro- 
cedure gives the orthonormal frame on R 3 

dx 2 dz' dy 

which generates the distribution 

V 2 

A = kcr(dz - ycte). 

Proposition 2. Consider the almost- Riemannian metric g = dy 2 + \dz 2 on 
R 2 . The y-axis is the union of the two geodesies starting at the origin with initial 
covectors (p y — ±l,p z — 0). The geodesies with initial covector (p y = ±l,p z = 
0) are given by 

/2 

V(t) = -Py^j=cn(K + t^\X\), 



m = "^^[tVW\+2sn(K + t^\)cn(K + t^\)dn(K + t^\X\)}, 

where K is the complete elliptic integral of the first kind 

r*' 2 dip 



S' 7 ' 

Jo 



1/2 sin 2 (p 



and cn,sn,dn denote the Jacobi elliptic functions of modulus k = Moreover 
the following properties hold true. 

(1) The almost- Riemannian spheres centered at the origin are subanalytic. 

(2) For A / 0, the cut point coincides with the first return to the z-axis that 
occurs att = 2Kj \f\\\, where two extremals with the same length intersect. 
The cut locus from the origin is {(y, z) \ y = 0} \ {(0, 0)}. 

(3) For A ^ 0, the conjugate point corresponds to t ~ ZK / ^f\\\. The conjugate 
locus from the origin accumulates at the origin as a set of the form {(y, z) \ 
z = ay 3 } U {(y, z)\z = -ay 3 } \ {(0, 0)} 7 with a £ 0. 

Proof. Define F 3 — and Pi =< p,Fi(q) >, i — 1,2,3. Using the Pontryagin 
Maximum Principle, the equations for normal extremals of the sub-Riemannian 
structure are given by the Hamiltonian system associated with H = \{P 2 + P 2 ), 
i.e. 

± = Pi, Pi = yP 2 P 3 , 
y = P 2 , P2 = -yPiPs, 
i = £p u P 3 = 0. 

There are three first integrals, namely p x ,Pz = \,H. The normalization condition 
H = 1/2 at t = gives 

Pi(0) 2 + P 2 (0) 2 = l, 

hence, we set Pi(0) = simp, ^2(0) = cost^. The set of extremals is invariant under 
the action of the group generated by the diffcomorphisms (x,y,z) 1— ^ (x, — y, z) 
and (a;, y, z) ^ (—x, y, —z). Therefore, it is sufficient to integrate the system with 
initial point (0, 0, 0) and covector (sin ip, cos ip, A) with A > and cos tp > 0. Recall 
that the normal extremals for the almost-Ricmannian structures are projections on 
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the (y, z) coordinates of the geodesies for the sub-Riemannian metric satisfying the 
transversality condition p x — 0. Hence, for A = we get y(t) = cos(ip)t, z(t) = 0. 
Assume now A > and set k, k' such that k 2 = 1 ~ s * nip , < k, k' < 1 and k 2 + k' 2 
1 . Then we find 

2 2 

y 2 = (l-P 1 )(l + P 1 ) = (l- Px -lLp z )(l + p x + lLp z ) 
= (2k 2 -^\)(2k' 2 + ^\). 
Setting y(f) = ^§y(t), the evolution equation for y is 

^(l-y 2 )(fc' 2 + fc 2 y 2 ), 
that can be integrated, with y(0) > 0, as y(t) = —cn(K(k) + t\~\, k), where 

«' 2 dtp 



K{k) 
Hence 



o \/l — k 2 sin 2 if 



2k 

y{t) = pen (K{k) + tV\, k). 

v A 

Remark that the extremals that project on geodesies for the ARS satisfy the 
transversality condition p x = which implies k 2 = 1/2. Thus the y coordinate 
of the geodesic with p x = Q,p y = 1, A > is 

y (t) = -^L cn (K + tVx), 

where, to simplify notations, we denote K(\/2/2) by K and we omit the dependence 
of the Jacobi function cn on the modulus. To compute the z coordinate along the 
same geodesic, we use the primitive J cn 4 (K + u) du = | [\u + sn (K + u) cn (K + 
u) dn (K + u)] (where k = \f2/2). Then 

z (t) = -^T^[tV\ + 2 sn(K + tV\) cn{K + tV\) dn (K + tVX)]. 

Using the symmetries of the system we find the required expressions for the geodesies 
starting at the origin for the almost-Riemannian metric. In particular, y and z are 
quasi-homogeneous with respective weights 1 and 3. The cut instant of a geodesic 
coincides with the first return to y = that occurs at t = 2K/VX, thus the cut 
locus is the z-axis. The conjugate time satisfies t ~ 3K / \/A, whence the conjugate 
locus can be approximated by the parametric curve 

V = 



A l/2> A 3/2 

(for the detailed proof see [1] ) . ■ 
Remark that for the desingularized structure, i.e., the sub-Riemannian Martinet 
flat case, the sub-analyticity of the sphere is lost in the abnormal direction for which 
k — > 1. This does not arise for the almost-Riemannian structure, since geodesies 
satisfy k 2 = 1/2. 

Property Q] of proposition [5] can be generalized to the generic tangential case 
using [11], see also the computations in section [5] 
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4.2. Analysis of the generic model of order 0. The objective of this section 
is to lift the generic model of order into a constant rank sub-Riemannian model 
in order to analyse the role of the invariants in the optimal dynamics. A geometric 
interpretation will be given in section [6] in terms of the form of the cut locus. 

Recall that from [7] the sub-Riemannian Martinet model of order zero is nor- 
malized to 

(1 + ayfdx 2 + (1 + fix + "iy) 2 dy 2 , 
where the distribution has the standard Martinet form 

2dz = y 2 dx. 

In this normal form the parameters a, (3, 7 are related to the geometric properties 
of the sphere with small radius and appear in the pendulum interpretation of the 
extremals. More precisely, for (3 — the extremal system is integrable while if /3 is 
non zero we have dissipation. In the integrable case the important parameter is a 
and if it is non zero the abnormal direction is strict. The role of the parameter 7 
is unimportant and it can be absorbed by reparameterization. 

Consider the almost-Riemannian metric on R 2 given by the orthonormal frame 

(42) Fl = ( £ z + y 2 /2 + £ y)^,F 2 = |-, 

where e 7^ 0. This metric can be seen as the generic model of order for an ARS 
in a neighborhood of a tangency point (use the normal form (F3) and weights of 
coordinates). Next proposition gives a possible lifting of the model of order at 
a tangency point for an almost-Riemannian metric, showing the relation with the 
model of order of a Martinet type distribution for a sub-Riemannian metric. 

Proposition 3. The generic model of order for an ARS in a neighborhood 
of a tangency point lifts into the sub-Riemannian Martinet model of order zero 
, ,f x yvi + 1 „ f v , — 7 — ttt on the distribution 2dz — y 2 dx = 0. 

(e(l+x))- ! (l+2e' y+o(y)) i « 

Proof. Applying the desingularization procedure (see section 13.2. 21) to the 
almost-Riemannian metric defined by (|42j) . we get the sub-Riemannian metric in 
R 3 defined by the orthonormal frame, still denoted by Fi,F 2 , 



dx 2 dz dy 



One gets 



[F 1 ,F 2 ] = -y(l + 3e'y) — , [[F 1 ,F 2 ] 1 F 2 ] = (1 + 6^) — , 

[[F 1 ,F 2 },F 1 } = sy(l + 3e'y)-^, 

hence the Martinet surface is the set {(x, y, z) \ y = 0}. Moreover the singular 
control in the Martinet surface is defined by 

ui det(F u F 2 , [[Fi,F 2 ],Fi]) + u 2 det(Fi, F 2 , [[F u F 2 ], F 2 ]) = 

which implies u 2 — 0. The corresponding trajectories are solutions of 

x = ui, y = 0,z — uiez. 

In order to build a coordinate system (x, y, z) in which the distribution has the 
normal form D = ker u, us = dz — y 2 /2dx, we normalize the singular flow to 
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lines parallel to the i-axis and lying in the Martinet surface. We consider the 
diffcomorphism 



x 




—e 



V = vV 1 + 2e'?y = y + e'y 2 + o(y 2 ), 
z = ze~ ex . 



The orthonormal frame in the new coordinate system becomes 

F 1 = -e(l + x) J, - e(l + 5)^ A, F 2 = (1 + 2e'y + (y)) A 



Hence the distribution is in the normal form dz = \dx and the metric is given by 

_ dx 2 dy 2 



■ v e 2 (l + £) 2 (l + 2 £ 'y + (y)) 2 ' 

■ 

Introducing F 3 = ^ and Pi ~< p,Fi >, the extremal flow is given by 

X = -e(l + X)P u 

Y = {l + 2e'Y + o(Y))P 2 , 

Z = -s(l + X)^-P 1 , 

Pi = -e(l + X)Y(l + 2e'Y + o{Y))P 2 P 3 , 

P 2 = e(l + X)Y{l + 2e'Y + o{Y))P 1 P 3 , 

P3 = 0. 

Setting P 3 = A and using the time parameter r such that dr = (1 + X)dt, we can 
write 

dX 

e?y (1 + 2e'Y + o(Y)) 



dr l + X 

dZ Y 2 



2, 



dPi 

~dT 

dP 2 
dr 



-XeY(l + 2e'Y + o(Y))P 2 , 
XeY(l + 2e'Y + o(Y))P 1 . 



Define 9 in R/27rZ by Pi — cos(6) and P 2 = sin(0). It satisfies 

^ = XeY(l + 2e'Y + o{Y)) 

and then 

d 2 i + 4 £ 'r + (y) . 
^ = Ae TTx sm(0) ' 

which can be approximated by 

^ = Ae(l - X + Ae'Y + o{Y)) sin(0). 
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According to [7] , this corresponds to a dissipative pendulum, the non nullity of the 
parameter e' inducing a coupling with the y-coordinate. Note that more computa- 
tions are necessary to get the sub-Riemannian Martinet metric in the normal form 
of order 0, leading to a true dissipative pendulum equation with no coupling with 
the x and y variables, see [7]. 

5. ASYMPTOTICS OF THE WAVE FRONT 

In this section we use the techniques and results from [7] , developped in the sub- 
Riemannian Martinet case, to compute asymptotics of the front from the tangency 
point for the generic model of order for ARS. Remark that the higher order 
terms in the expansion of the elements of the orthonormal frame play no role in the 
estimation of the front and, consequently, in the estimation of the cut and conjugate 
loci from the tangency point (see section [BJ, as one can check easily. 

Proposition 4. Consider the ARS on R 2 defined by the orthonormal frame 
given in \4->ty - The extremals satisfying initial condition 

{y, z,p y ,Pz)\t=o = (0, 0, ±1, A) 

with |A| ~ +oo can be expanded as 

y(t) = rjY^t/rti + rfY^t/rfi+otf), 

z{t) = rfZ {tlri) + ^Z\t/ri) + o{rf), 

where rj = 



(53) y° = p°, p° = - [ z) y ' 



with initial condition (Y°, Z°,P Y , P|)|t=o = (0,0,±1,±1) and 
(54) Y 1 = Pi, 

Z 1 = ^PUY^f + P^iYyY'+eZ^Y^+e'iY ) 5 ), 

P X y = -PlPh(Yy-{P^ z )\\{Y^+eZ^ + \e'{Y^)% 

Pi = -\(P°z) 2 e(Y°) 2 , 
with initial condition (Y 1 , Z 1 , Py , P^)|t=o — (0,0,0,0). 

Remark. Computations in this case are similar to the ones of the Martinet sub- 
Riemannian case. System (I54[) represents a variational equation whose integration 
is related to the second-order equation 

Y l + {-PfY° 2 )Y l = K(Y°) 

where Y° is a periodic elliptic function. 

Proof. Recall that y, z have weight 1 and 3, respectively. In order to have the 
standard Darboux form of order 1, we fix the weight for p y and —2 for p z . 

The Hamiltonian is 

H =\{pl{ez + V y2 + e'y i f+pl) 



and the extremal flow is 

V = Py, Py = -pl{£z + y 2 /2 + e'y z ){y + ?>e'y 2 ), 

z = p z (ez + y 2 /2 + e'y 3 ) 2 , p z = -p 2 z {ez + y 2 /2 + e'y 3 )e. 



SPHERE AND CUT LOCUS AT TANGENCY POINTS FOR 2-D ARS 



11 



Y 



According to the weights we set 

V = V Y , Py = P Y, 
z = rfZ, p z = |f-, 

where r\ is a parameter. The evolution equations for (Y,Z,Py,Pz) are 

Py 
V 

Y 4 

Z = P z (— +eZY 2 + e'Y 5 +e'TiY e + 2ee'-nZY 3 + £ 2 r]Z 2 ) 7 

Py = -Pl(^ +eZY +^-e'Y 4 + 3r 1 e'Y 2 (eZ + s'Y 3 )), 
Zrj 2 

Pz = -Ph^+evZ + e'rjY 3 ). 
Considering the expansions with respect to 77 

Y = Y° + TJY 1 +0(77), P Y = P y + r]P Y + 0(77), 
Z = Z + V Z 1 +o( V ), P z =P z +r)P z +o( V ), 

by identification we find that the leading terms satisfy 

y = Pi p = (P%f(Y°f 

V ' , r Y 2»7 ' 

z ° = p °z = 0. 

In particular P z is constant. Setting A = p z (0), for A ^ we can fix 77 = l/^/[X[ 
and then P z is normalized to 1 or -1. 

Introducing the time parameter s = iy^Aj the equations ()55[) for the first-order 
approximation become 

dY°_ _ p0 dP°_ _ (Pg) 2 (r°) 3 

ds ~ r Y> ds ~ 2 ' 

dZl = P°,(Y f p a ±1 
ds 4 ' Z 

System (|56[) coincides with the Hamiltonian system for the nilpotent model that 
has been integrated in Proposition [21 using elliptic functions with modulus k such 
that k 2 = 1/2. The solution is given in Proposition^ 

Y°(s) = -P Y (0)V2cn(K + s), 
P° 

Z°(s) = -^-(s + 2six(K + s)cn(K + s)dn(K + s)), 
P Y (s) = Py(Q) + (P Y (0)) 3 (-l + V2dn(K + s)sn(K + S )), 
*°z( 

Using s — t/r), the system for (Y, Z,Py,P z ) becomes 
dY Py, 



P° z (s) = ±1. 



ds 
dZ 

Is '""A 



,,Z Pz (Xl + V ( £ ZY 2 + e'Y 5 ) + v 2 (e' 2 Y 6 + 2ee' ZY 3 + e 2 Z 2 )), 



dPs 



- = -P 2 z (—+fi{eZY+yY 4 ) + ii 1 2 e'Y 2 {eZ + e'Y 3 )), 



3 

ds ' z ^~2 ' ' 2 

^ = -P 2 ze(r^+ V 2 {eZ + e'Y 3 )). 
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Hence, identifying terms of order 0, one gets 



=- = ipifr°) 4 + i«((y ) 3 y 1 + ez (y ) 2 + e , (y ) 8 ), 

as 4 

^ = -P z pi(Yy-(P z )^(Y ) 2 Yi+eZ°Y + y(Y )% 
ds 2 z 



6. Geometric estimates of the conjugate and cut loci 

6.1. The conjugate locus. The following result gives a description of the conju- 
gate locus from a tangency point of a 2-dimensional ARS. 



Proposition 5. Consider an ARS on R defined by the orthonormal frame 

' dz ' dy 



F\ = (ez + y 2 /2 + e'y 3 + o 3 (y, *))■£-, F 2 



T/ien </iere exists a constant a^O swc/i f/iaf the conjugate locus from (0, 0) accu- 
mulates at (0, 0) as the set 

{(y, z)\z = ay 3 } U {(y, z) \ z = -ay 3 } \ {(0, 0)}. 

Proof. Applying Proposition 21 the exponential map at (0,0) is given by 

(? y , s) i ^ (t^s) + o{ V ),n 3 Z\s) + o( V 3 )), 



where s = ty / p z (0), 77 parametrizes the initial covector as (p y {0) = ±l,p z (0) = 
P°/t7 2 ), and 

Y°(s) = -P^(0)V2cn (K + s), 

Z°(s) = P z ±(s + 2sn (K + s)cn (if + s)dn (K + s)). 

The conjugate time is the first zero of the Jacobian of the exponential map. The 
Jacobian is equal, up to a multiplicative constant, to 

HZ r\V° 

r/ 3 (y od^_ 3Z ody 
ds ds 

It was proven in [7] that the function 

«•>-*•<•>£ -•* ,! ir 

has its first positive zero at s = so ~ 3K and that j'(so) 0. Hence, the conjugate 
time is of the form sq + o(l) where o(l) is a continuous map going to zero when 77 
goes to zero. 

In terms of the singularity theory, this computation proves that the exponential 
map of a general two-dimensional ARS can be seen as a small deformation of 
the exponential map of the nilpotent case. In the nilpotent case, the exponential 
map has only stable singularities (folds) corresponding to the first conjugate locus. 
Hence, in the general case, the first conjugate locus also corresponds to folds and 
accumulates at (0,0) as the set {(y, z) \ [z — ay 3 ){z + ay 3 ) — 0}, where q^O, 
see Proposition [2] ■ 
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6.2. The cut locus. In this section we provide a description of the cut locus at a 
tangency point for a generic ARS. As one can infer from the proof of the following 
proposition, the shape of the cut locus is determined only by the terms of order 
up to in the expansion of the elements of the orthonormal frame. Higher order 
terms do not contribute to the estimation of the way the cut locus approaches to 
the tangency point. One can see in figure [5] small spheres for different values of e 
and e' . 




Figure 2. The sphere of small radius for the nilpotent approxi- 
mation (dotted line) and for an example with e' = (dashed line) 
are symmetric. The two spheres are not C 1 at their intersection 
with the cut locus, which in both cases is the vertical axis. The 
sphere of small radius for the generic model of order zero in which 
e' (solid line) loses the symmetry. In this case, the cut locus 
is different from the previous cases (see Proposition 6). 



Proposition 6. Consider the ARS on R 2 defined by the orthonormal frame 

Fi - (ez + y 2 /2 + e'y 3 + o 3 (y : z)) — ,F 2 = — . 

oz ay 

Then, if e 1 ^ 0, there exist non zero constants a.\,a.i such that the cut locus from 
(0,0) accumulates at (0,0) as the set 

{(y, z)\z>0,z 2 - a iy 3 = 0} U {(y, z)\z<0,z 2 ~ a 2 y 3 = 0}. 

Proof. In the following, we restrict our analysis to the upper half plane z > 0, 
the computations being equivalent in the case z < 0. First of all, recall that for the 
nilpotent model (e = e' = 0), the cut locus is {(y, z) | y = 0} \ {(0, 0)} and the cut 
time on the geodesic with initial covector (1, A) is 2K/y r X which corresponds to the 
first intersection with the symmetric geodesic whose initial covector is (—1, A), see 
Proposition O 
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Denote by (Y°, Z°, P Y ,P Z ) the geodesic with initial condition (0, 0, 1, 1) for the 
ARS with e = e' = 0, i.e., 

Y°(s) = -V2cn(K + s), 

Z°{s) = -(s + 2sn(K + s) cn(K + s) dn{K + s)), 

Py(s) = \/2 dn (K + s) sn (K + s), 
Pi = 1. 

Moreover, denote by Y X ,Z ,P Y ,P Z the terms of order in the expansion of the 
geodesic for the ARS with e' — 0, i.e., solutions of the system ([54"f with e' = and 
initial condition (0, 0, 0, 0). Set (y 1 , Z 1 ,Vy , Vz ) to be the terms of order in the 
expansion of the geodesic for the ARS with e' 7^ 0, i.e., solutions of system (j54")l 
with the initial condition (0,0,0,0). Finally, define four functions of s, 31, 52, 33 
and 34, by 

y 1 =Y 1 +e f g 1 , Ty 1 = Py 1 + e 'g 3 , 
Z 1 = Z 1 +e , g 2 , Vz 1 = Pz 1 +e'g 4 . 

Combining the equations satisfied by Y 1 , Z 1 , P Y , P\ and by y 1 , Z 1 , Vy 1 . Vz 1 , 
we find that (31,32,33,34) satisfy the following system of ODEs 

91=9* 93 = -l(Y°) 2 9i-l(Yy 

32 = 9i(Y ) 3 + (Y ) 5 34^0, 

where 3^ = dgi/ds, with the initial conditions 31 (0) = 32(0) = 33(0) = 0. 
Remark that if 

(Y°,Z°,P Y ,P z ,Y 1 ,Z\P Y ,P 1 z ,g 1 ,g 2 ,g 3 ) 

is solution of f(p ]) .([M | .([67 ]l ) with the initial condition (0,0,1,1,0,0,0,0,0,0,0) 
then 

(-y°, z°, -p Y , p° z , -y\z\-p y , p z , gi , -32,33) 

is also solution with the initial condition (0, 0, —1, 1, 0, 0, 0, 0, 0, 0, 0). Moreover, one 
can compute numerically that g\(2K) ~ — 2tt, g2(2K) ~ — tt and 33(2A') ~ 0. 

Let us compute the front at time t = 2Krj close to the initial condition rj , that 
is for t] = rjo + ct/q + 0(77^). Making Taylor expansions in terms of 770, one finds for 
the front corresponding to the initial conditions p y (0) = 1 and p z (0) = 1/ rj 2 

V = r 1 Y (^) + V l{Y^2K) + e'g l {2K)) + o{ V l) 1 
z = ^Z°C-^o)+ V ^(2K) + e'g 2 (2K))+o(^). 

Hence we obtain 

y = gl{Y\2K)+e>g 1 {2K) + 2Kc) + o{r 1 l) 1 

z = V $Z°(2K)+r 1 l(Z 1 (2K)+e'g 2 {2K) + 3cZ Q (2K))+o{4) : 

since Y°{2K) = Z°(2K) = and Y°{2K) = -1. For the front corresponding to 
the initial conditions p v (0) = — 1 and p z — 1/fj 2 where fj = r/o + c'rj 2 + o(g 2 ) one 
finds 

y = V 2 (-Y 1 (2K)+e'g 1 (2K)-2Kc') + o( V 2 ), 

z = T 1 IZ°(2K) + ^(Z 1 (2K)~e'g 2 {2K) + 3c'Z°(2K))+o(4). 

These expressions are afline with respect to parameters c and c', up to order 2 
for y and 4 for z in the variable t]q. The two geodesies with initial covectors 
p y (0) = l,p z (0) = I/77 2 andp a (0) = -l,p z (0) = 1/fj 2 intersect for 

. / _ Y 1 (2K) 



-R ro(l), 

d-c = e ' 92 £ K) + o(l) 
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where o(l) denotes any function going to with 770- Hence the intersection is for 
and 

, e'g 2 (2K)-Y 1 (2K) 

c= fie + 

which implies that the cut point is 

ycut = v 2 oe'(9i(2K)-g 2 (2K)) + o(r 1 i), 
Zcut = Vo^f + o(vo)- 

Hence, if e' 7^ 0, the upper branch of the cut locus from (0, 0) accumulates as the 
set {(y, z) I z > 0, z 2 = otiy 3 }, where 

AK 2 AK 2 



ai 



9e' 3 (5i(2if ) - g 2 {2K)f 9e' 3 7r 3 ' 



Similar computations show that the lower branch of the cut locus from (0, 0) accu- 
mulates as the set {(y, z) \ z < 0, z 2 — a 2 y 3 }, where 

AK 2 AK 2 

Q/rj r-O 

9£' 3 ( 5 i(2^)-|-(? 2 (2^)) 3 3V 3 tt 3 ' 

■ 

Remark that the case of generic ARS with e' 7^ is rather different from the 
sub-Ricmannian Martinet case (see [7] ) , in which a similar argument cannot apply. 
Indeed, in the latter situation, the asymptotic expansions for small time cannot 
be used, since there exists an abnormal direction corresponding to the case where 
k — > 1 and K (k) — > 00. Hence, we should use the asymptotic expansion for a time 
parameter tending to +00, which is clearly not valid. 
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